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We explore the effects of noncommuting applied fields on the ground-state ordering of the quasi- 
one-dimensional spin- 1/2 XY-like antiferromagnet CS2C0CI4 using single-crystal neutron diffraction. 
In zero field interchain couplings cause long-range order below Tm~217(5) mK with chains ordered 
antiferromagnetically along their length and moments confined to the {b, c) plane. Magnetic fields 
applied at an angle to the XY planes are found to initially stabilize the order by promoting a spin-flop 
phase with an increased perpendicular antiferromagnetic moment. In higher fields the antiferromag- 
netic order becomes unstable and a transition occurs to a phase with no long-range order in the 
(6, c) plane, proposed to be a spin liquid phase that arises when the quantum fiuctuations induced 
by the noncommuting field become strong enough to overcome ordering tendencies. Magnetization 
measurements confirm that saturation occurs at much higher fields and that the proposed spin-liquid 
state exists in the region 2.10< Hsl <2.52 T || a. The observed phase diagram is discussed in terms 
of known results on XY-like chains in coexisting longitudinal and transverse fields. 

PACS numbers: PACS numbers: 75.25.+Z, 75.10.Jm, 75.30.Cr, 75.45.+j 



I. INTRODUCTION 

A promising area for the study of zero-temperature 
quantum phase transitions is that of quantum spin sys- 
tems in noncommuting applied magnetic fieldsa. Such 
noncommuting terms introduce quantum fluctuations 
into the T — ground state which, for large enough 
fields, can completely disorder the system. Such a situa- 
tion has been studied in considerable detail in the three- 
dimensional (3D) laing ferromagnet LiHoF4 in a trans- 
verse magnetic fieldo. However, due to its high dimen- 
sionality, the system behaves in a mean-field like way. In 
this paper we consider a one-dimensional (ID) quantum 
magnet and the effects of a noncommuting field on its 
ground state. 

The scenario that we investigate experimentally is that 
of the ID spin-1/2 XXZ model {Hxxzi^)) in a noncom- 
muting field (Happ), given by the Hamiltonian, 

n = nxxz{A) + nAPP, (i) 

where, 

HxxziA) = J ^ (5f 5f+i + SfSf^, + A5f 5f+0 

i 

(2) 

and 

nAPP^Y.^'^'^ + ^^^^^ (3) 

i 

and where J > is the antiferromagnetic (AF) ex- 
change constant, and A < 1 is the anisotropy param- 



eter (a gx,zl^B factor is incorporated in the magnetic 
field B^-^). The operators S", where a — x,y, z, are the 
usual spin operators for spin-1/2, and the applied field 
term does not commute with the exchange Hamiltonian 
[nxxz{A),nAPp] ^ for |A| ^ 1 and ^ 0. Indeed, 
the system we shall study, CS2C0CI4, is an excellent re- 
alization of Ti.xxz{A) where A = 0.25, and therefore 
should Ue well approximated by the famous XY model 
(A = 0)0. 

The physics of the TCxxziA) model at A = is well 
known in the absence of a magnetic field i.e. = B^ = 
0. It is that of a ID noninteracting Fermi gas as shown 
by the Jordan- Wigner transformationQ. The effect of a 
field along z, B^ ^ 0, (commuting) is trivial; it acts 
as a chemical potential and changes the filling level of 
the fermions in the chain inducing magnetization along 
z. Below a critical field, J5^, representing the field of 
complete saturation of the chain, the excitation spectrum 
remains gapless and the T = correlation functions fall 
algebraically as power laws. In contrast, the physics for 
the case of a field along B^ ^ 0, B^- = 0, is neither 
trivial nor widely known. 

The action of noncommuting field B^ on 7ix^ft(0) has 
been considered theoretically by Kurmann et a/.u'u. From 
these studies it was found that the in-plane field has two 
effects: 1) it breaks the U(l) symmetry of the XY-model 
to a lower, Ising-like, symmetry which causes the ground 
state to long-range order (LRO) at T = into a spin-flop 
type Neel state. In fact at a special coupling, B^ = \f2J 
the spin-flop Neel state is the exact T = ground state. 
2) The second effect is to introduce quantum fluctuations 
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FIG. 1: Scliematic representation of tlie pliases of tlie ground 
state of the XY-model as a funejijipn of applied transverse field 
as proposed by Kurmann et al.WB (see text for details). In the 
absence of a field the correlation functions decay as a power 
law. Small magnetic fields induce perpendicular long-range 
ordered antiferromagnetism and the ordering is characterized 
as a spin fiop phase. At a field \/2J the classical Neel state 
is the exact T = ground state. The antiferromagnetism 
is rapidly suppressed at higher fields by the quantum fiuc- 
tuations induced by the noncommuting field. Above about 
1.5J long range order is destroyed by these fiuctuations and 
the ground state is characterized as a spin liquid (S.L.) state 
with exponentially decaying correlations in the spin compo- 
nents perpendicular to the field. Above a crossover field of 
about 2 J nearly all the spin moments are aligned along the 
field direction and the physics here is characterized as a sat- 
urated phase. 



that of magnetic chains which are ordered antiferromag- 
netically along their lengths with ordered moments lying 
in the (6,c) plane. The experimental results show that an 
applied field initially further stabilizes the ordered struc- 
ture (spin- flop phase), but as the field increases fluctu- 
ations induced by field noncommutation cause a sharp 
transition to a new phase which we propose to be the 
spin-liquid phase discussed above. Study of the ferro- 
magnetic component ( III DD confirms that the saturation 
field is higher still, and that the proposed spin liquid 
phase exists in the region 2.10 < Hsl < 2.52 Tesla \\a. 

An extended discussion of our results is given in 
Section IV. The observation of commensurate order 



(=(0,1/2,1/2)) in CS2C0CI4 i s rela ted to the effective 
Ising-type interchain coupling ( IV A ) and coropared with 
the quasi-elastic results of Yoshizawa et al xi A micro- 
scopic study of the ordering is presented in |IV B| based 
on a mean-field analysis of the ground-state energy. The 
mag nitude of the ordered mom ent is established in Sec- 
tion |IVC| . Finally, in Section |IVD| the magnetization 
curve is critically examined and compared with known 
results for anisotropic spin chains in noncommuting ap- 
plied fields. 



II. EXPERIMENTAL DETAILS 



A. Properties of CS2C0CI4 



into the system. At high fields this causes a phase transi- 
tion to occur where the fluctuations become large enough 
to destroy the LRO altogether. This disordering field is 
below that where the system reaches its saturation point. 
This phase transition is therefore a nontrivial quantum 
phase transition through a quantum critical point with 
the noncommuting field as a control parameter. Fig. |l| 
shows a schematic outline of the physics of the XY-model 
in a transverse field. j.. j.. 

The theoretical studiea3'El suggest that this behavior 
is generic to a wide class of magnets in noncommuting 
fields, and a disordered spin-liquid (S.L.) phase is ex- 
pected for all Hxxz{\'^\ 7^ 1) in equation (1). Spin Hq- 
uid phases such as this are generally gapped (evidenced 
by the exponentially decaying correlations in the zero- 
temperature ground state) and therefore robust against 
small perturbations. We therefore expect the effects of 
nonzero A and interchain coupling to modify the tran- 
sition fields in Fig. |l| but not to change its qualitative 
content. 

In this paper we present a detailed study of the ground 
state ordering of the quasi-lD XY-like antiferromagnet 
CS2C0CI4 in noncommuting fields from zero to well above 
saturation. The crystalline and magnetic properties of 
CS2C0CI4 and the experimental tool of neutron diffrac- 
tion will be introduced in Section |l[ 

Section III presents the observed commensurate mag- 
netic structure and its field dependence. The structure is 



CS2C0CI4 has been proposed as a spin-1/2 quasi- 
lD XY-like antiferromagnet with chains running along 
the 6-direction. Heat-capacity measurementsH showed a 
broad maximum (characteristic of low-dimensional sys- 
tems) around T ~ 0.9 K and the overall temperature 
dependence agreed well with numerical predictionsLl for 
a 5=1/2 XY-like AF chain with exchange coupling J = 
0.23(1) meV. A small lambda-type anomaly observed 
in the specific-heat at T/y = 222 mK was interpretedEl 
as indicating a phase transition to a magnetically or- 
dered phase caused by small coupli ng s b etween chains. 
The in-plane magnetic susceptibilityliO'til also indicated 
magnetic ordering below 222 mK and the temperature 
dependence between 40 mK and 4.2 K was in excel- 
lent agreement with numerical calculations for AF XXZ 
chains with A = 0.25 in Eq ^ _ 

Quasi-elastic neutron scattering experimentsQ showed 
that between T = 0.3 and 0.6 K the critical scattering of 
CS2C0CI4 is sheet-like perpendicular to the b axis, thus 
confirming the proposed quasi- ID character of the mag- 
netic properties. Interestingly, the critical scattering is 
driven partly to incommensurate positions due to com- 
peting interchain interactions. Those earlier diffraction 
experimental concentrated on the critical scattering at 
temperatures above the ordering transition of 0.22 K in- 
ferred by macroscopic measurements. As the ordering in 
this material promises to have some unusual and chal- 
lenging features, we undertook detailed neutron diffrac- 



3 




FIG. 2: Crystal structure of CS2C0CI4. The figure shows 
12 Co'^^ ions, each surrounded by a distorted tetrahedron of 
Cl~-ions. The Co^^ ions interact mainly via an AF super- 
exchange interaction J along the b axis, forming AF spin 
chains which interact weakly via super-exchange interactions 
Jab, Jac and Jbc as explained in the text. The shaded rect- 
angular planes indicate possible orientations of the XY easy 
planes following Ref. hsl 



tion experiments to determine the magnetic structure 
and the behavior in applied noncommuting magnetic 
fields at temperatures much below the proposed order- 
ing transition of 0.22 K. 

CS2C0CI4 crystallizes in the orthorhombic (apd non- 
symmorphic) space group Pnma {D\^^ No. 62)EHl. The 
crystal structure is shown scbjaaatically in Fig. |[ The 
lattice parameters at 0.3 K ardj-B a = 9.71 k ,b = 7.27 k 
and c = 12.73 A. The magnetic ions, Co^+ wittpSpin 
5 = |, occupy site 4c in the unit cell at positional^llJ 

(1) : (0.235,0.25,0.422) 

(2) : (0.735,0.25,0.078) 

(3) : (0.765,0.75,0.578) 

(4) : (0.265,0.75,0.922) . 

Neighboring spins interact via a superexchange interac- 
tion involving a bridge of two Cl~-ions with the path 
Co2+-Cr-Cr-Co2+. The shortest Cr-Cl- distance is 
between neighbors along the 6-axis separated by 3.63 A. 
As this is close to twice the ionic radius of CI a sizeable 
overlap of electron wave- functions contributing to the ex- 
change integral is expected. 

Among other possible exchange paths are Jac between 
sites 1 and 2 (dci--ci-= 4.05 A), Jab between sites 1 
and 3 (rfci^-ci"^ 4.04 A), and Jhc between sites 1 and 
4 =4.01 A). Since the overlap of the electronic 

wave-functions decreases very rapidly (approximately ex- 
ponentially) with distance, these exchange paths are ex- 
pected to yield much smaller interactions than the cou- 



pling J along b such that CS2C0CI4 can be regarded as a 
system of weakly-coupled spin chains along the h axis. 

Each Co^+ ion is tetrahedrally coordinated by Cl^ lig- 
ands. Small distortions from a perfect tetrahedron lead 
to a splitting of the S—^ orbital singlet state into two 
Kramers doublets with a separation 2£'=1.3(1) meV. 
The magnetic exchange energy is much lower than the 
inter-doublet separation and therefore only the lowest- 
lying doublet states participate in the low-energy dynam- 
ics at low temperatures (T <C I?). Projecting the Heisen- 
berg exchange between the true spins onto the lowest- 
lying doublet of I ± i) states gives an effective spin-1/2 
Hamiltonian with XY-like exchange, Tixxz in Eq- § with 
A=0.25. 

Rotations of the C0CI4 tetrahedra in the unit cell lead 
to different orientations of the XY easy plane between 
sites (1,3) and (2,4) and give the 6-axis as the only com- 
mon in-plane direction for all sites. The main distortion 
of the C0CI4 unit from a perfect tetrahedron is due to 
one of the four Cl^ ions being rotated by several de- 
grees around the 6-axis with respect to the central 00^+ 
ion and previous studiescJ proposed that the normal to 
the XY easy plane (z-axis) bisects this large angle. This 
gives the z-axis direction (sin/3,0,cos/3) on sites (1,3) and 
(-sin/3,0,cos/3) on sites (2,4) with /3 = -38.8°. The 
shaded rectangular planes in Fig. ^ show the XY planes 
in this case. Another possibility is that the z axis is along 
the vector connecting the central Co^+ ion with the Cl^ 
ion rotated most and in that case /3=+19.4°. 



B. Experimental method 

A high-quality 6.45 g single crystal of CS2C0CI4 was 
grown from solution. The crystal was aligned with its 
(0, fc, I) plane in the horizontal scattering plane and was 
cooled to temperatures between T = 80 and 250 mK 
using an Oxford Instruments dilution refrigerator insert 
placed inside a vertical 7 Tesla superconducting magnet. 

Neutron diffraction measurements were made at the 
Hahn-Meitner Institut in Berlin, Germany. The two-axis 
crystal diffractometer E6 was employed with an incident 
energy of Ei — 14.72 meV and with a pyrolytic graphite 
(PG) monochromator in double-focusing mode to in- 
crease the neutron flux at the sample position. The scat- 
tered neutrons were counted in a 20° wide BF3 detector 
bank with position sensitivity along the horizontal direc- 
tion giving 200 channels. This gave an angular resolution 
of 0.1° in the total scattering angle 28. The intensities 
of nuclear and magnetic Bragg reflections were measured 
as a function of both 20 and the sample rotation an- 
gle ^I^, thus constructing full two-dimensional (2D) maps 
of the scattering intensity in the (5*, 20) plane. This al- 
lowed simultaneous coverage of both magnetic signal and 
background. Typical counting times were 5 minutes per 
2D map to determine the total integrated intensity, and 
40 seconds for a 20 scan at the peak center to extract 
the peak intensity. A two-dimensional Gaussian with ad- 
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FIG. 3: Integrated intensity of the antiferromagnetic 
(0, 0.5, —1.5) reflection vs. temperature in zero external field. 
Intensity units are the same as in Table ^ The solid line is a 
guide to the eye. The inset shows the Bragg peak intensity as 
a function of ^ (left) and 2Q (right) at temperatures below 
(solid circles) and above (open circles) the transition temper- 
ature Tn~217 mK. Solid lines come from a two-dimensional 
fit to the data in the {^,20) plane as described in the text. 



justable rotation of the main axes of the ellipsoid gave a 
good account of the observed peak line shapes in both 
and 28 and the intensities of the Bragg reflections were 
obtained from least-square fits to the experimental data 
to reduce the sum of discrepancies x^. 

III. EXPERIMENTAL RESULTS 
A. Magnetic order in zero field 

Upon cooling below T/v=217 mK extra Bragg reflec- 
tions were observed at the commensurate (0, n -|- 0.5, to -|- 
0.5) reciprocal lattice positions with n and m integers, in- 
dicating a transition to a magnetically long-range ordered 
state. Fig. |3| shows the temperature dependence of the 
(0, 0.5, —1.5) AF reflection showing the onset of order be- 
low Tjv. The extracted transition temperature Tn=2\7 
mK is consistent with that infeped from specific-heat and 
susceptibility measurements. eIEj 

The observed magnetic reflections are associated with 
a magnetic ordering wavevector k — (0, 1/2, 1/2). To de- 
termine the magnetic structure we flrst used group the- 
ory to identify the spin configurations consistent with the 
wavevector k for the given crystal symmetry, and second 
we compared the structure factor of possible spin configu- 
rations with the experimentally observed magnetic Bragg 
peak intensities. 

The group theory analysis and determination of the 
symmetry-allowed basis vectors is presented in Appendix 
A. After comparison of the data with the possible eigen- 
vectors, we find that the observed structure belongs to 
the r^" irreducible representation with eigenvector 0^" 
given in Eq. nfl. This eigenvector has 6 degrees of free- 
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FIG. 4: Magnetic structure of CS2C0CI4. Spins (indicated by 
arrows) order antiferromagnetically along chains (shown by 
dashed lines in (a)). Ordered moments are contained in the 
(6, c) plane and make a small angle with the fe-axis (see text for 
details). Relative ordering of the chains leads to degenerate 
domains (a)-(d) belonging to the same irreducible representa- 
tion with eigenvectors given in Eq. |l^. Labels 1-16 in (a) 
indicate the 16 spins in the magnetic unit cell (1-4 label the 
four atoms in the chemical unit cell shown in Fig. ^). Solid 
and open circles are -ions with height along the a-axis 
close to 0.25 and 0.75, respectively. 



dom corresponding to the three components of the mo- 
ments mi and TO3 . Using spherical coordinates these can 
be written as toi = Mi(sint?i,cost?i cos cos t^i sin(/)i) 
and TO3 = M3(sini?3,cosi?3cos03,cosi?3sin03). For 
"^x.z — spins are in the (6, c) plane and is the az- 
imuthal angle with the 6-axis. 

Fig. |4|(a) shows a pictorial representation of the 
eigenvector in Eq. |l^ in the special case of ordered spins 
contained in the (6, c) plane, making a small angle with 
the 6-axis (^i = — ^s) and having equal magnitude on 
all sites (|toi| — |to3|). The structure can be described 
in terms of antiferromagnetic chains along h with a cer- 
tain ordering pattern between adjacent chains. Starting 
with the basic structure shown in Fig. ^(a) other distinct 
domains shown in Fig. ^(b)-(d) can be constructed by 
changing the sign of either the h- or the c- spin compo- 
nents, or the relative phase between chains 1 and 3 (for 
details see Appendix A). 
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FIG. 5: Experimentally observed integrated magnetic Bragg 
peak intensities (open circles) fitted to the model for the mag- 
netic structure described in the text (solid line). Intensity 
units are the same as in Fig ^. The horizontal axis indicates 
the number of the Bragg peaks in the list given in Table ^ 
Inset: Schematic diagram of the {b, c) reciprocal space indi- 
cating the positions of magnetic Bragg peaks (solid circles) 
observed in the present experiment in the ordered phase be- 
low Tjv=217(5) mK. Grey circles indicate reciprocal lattice 
positions in the (6, c) plane and undulating solid lines mark 
the positions of the sheets of quasii-elastic magnetic scatter- 
ing observed in earlier experimentaJ at higher temperatures 
(r>0.3 K) above the transition to the ordered phase. 



We find that the measured Bragg intensities can be 
consistently described by an equal population of A and 
B domains of T^'^. The best fit of the model to the mag- 
netic intensities (for details see Appendix) is shown in 
Fig. ^ and in Table |, and gives a good description of the 
experimental data. A single-domain structure of either 
A- or B-type as defined in Appendix A and shown in 
Fig. H is unable to account for the results. 

Assuming spin moments confined to the (6, c) plane 
(i9i = i?3 = 0) and equal ordered moments on all sites 
(|mi| ~ \iTi3\) the best fit results are <j)i — 15(5)°, 
(/)3 = —15(5)° and fraction of A-domain a — 0.48(3) 
with the sum of discrepancies = 2.55. The obtained 
value for a ^ 1/2 demonstrates that domains A and B 
occu py th e sample in equal parts. It will be shown in Sec- 
tion pVB that those two domains have the same mean- 
field exchange energy and are thus expected to occur with 
the same probability. Taking mi out of the (6, c) plane by 
choosing a series of non-zero values for the out-of-plane 
angle di and fitting z?3 leads to worse agreement with 
the observed intensities, so within the accuracy of the 
experiment it was concluded that magnetic moments are 
mostly contained in the (&, c) plane. 

Fig. H compares the temperature dependence of the re- 
duced intensity ///T=80mK of three AF reflections. Here 
/T=80mK is the observed intensity at base temperature 
T — 80 mK. The three reflections have essentially the 
same temperature dependence, indicating that the mag- 
netic structure is unchanged between 80 mK and T^r. 



TABLE I: Experimentally observed magnetic Bragg peak in- 
tensities compared to the calculated intensities based on the 
model explained in the text. The magnetic Bragg reflections 
are in the order of decreasing observed intensity. A visual 
comparison between experiment and model is made in Fig. ^ 



Number 




Exp. Int. 


Calc. Int. 


1 


(0,0.5,1.5) 


55.1(1.8) 


51.4 


2 


(0,-0.5,-1.5) 


49.4(1.8) 


51.4 


3 


(0,0.5,-1.5) 


42.4(1.7) 


44.8 


4 


(0,-0.5,-4.5) 


16.2(1.5) 


15.1 


5 


(0,0.5,-4.5) 


15.5(1) 


15.5 


6 


(0,0.5,0.5) 


12.1(1.4) 


9.7 


7 


(0,1.5,1.5) 


10.4(0.8) 


10.1 


8 


(0,0.5,-0.5) 


10.3(0.6) 


8.1 


9 


(0,1.5,-1.5) 


10.1(0.9) 


8.4 


10 


(0,-1.5,-4.5) 


9.9(0.2) 


10.1 


11 


(0,1.5,-4.5) 


9.1(0.2) 


8.8 


12 


(0,-0.5,-3.5) 


7.8(0.2) 


8.8 


13 


(0,-0.5,-2.5) 


7.5(0.4) 


6.0 


14 


(0,0.5,-3.5) 


6.9(0.2) 


8.0 


15 


(0,0.5,-2.5) 


6.1(0.2) 


6.0 


16 


(0,0.5,-6.5) 


4.1(0.4) 


3.5 


17 


(0,2.5,-4.5) 


4.0(0.2) 


4.3 


18 


(0,1.5,-3.5) 


3.8(0.4) 


4.1 


19 


(0,-0.5,-6.5) 


3.3(0.4) 


3.4 


20 


(0,1.5,2.5) 


2.8(0.6) 


2.4 


21 


(0,2.5,-1.5) 


2.5(0.4) 


2.2 


22 


(0,3.5,-4.5) 


2.4(0.4) 


1.8 


23 


(0,2.5,3.5) 


2.2(0.6) 


1.9 


24 


(0,1.5,-0.5) 


2.0(0.6) 


0.84 


25 


(0,3.5,-3.5) 


2.0(0.6) 


0.61 


26 


(0,2.5,1.5) 


1.9(0.4) 


3.0 


27 


(0,1.5,-2.5) 


1.3(0.4) 


2.0 


28 


(0,2.5,-2.5) 


1.14(0.2) 


0.64 


29 


(0,-0.5,-5.5) 


1.14(0.4) 


1.2 


30 


(0,1.5,0.5) 


1.0(0.2) 


1.1 


31 


(0,2.5,2.5) 


0.95(0.4) 


0.8 


32 


(0,2.5,-3.5) 


0.94(0.2) 


1.5 


33 


(0,0.5,-5.5) 


0.84(0.3) 


1.1 


34 


(0,2.5,-2.5) 


0.79(0.6) 


0.64 


35 


(0,2.5,0.5) 


0(0.4) 


0.40 


36 


(0,3.5,0.5) 


0(0.4) 


0.18 


37 


(0,3.5,-1.5) 


0(0.4) 


0.83 


38 


(0,4.5,-0.5) 


0(0.4) 


0.08 


39 


(0,4.5,-1.5) 


0(0.4) 


0.39 


40 


(0,4.5,0.5) 


0(0.4) 


0.09 


41 


(0,4.5,1.5) 


0(0.4) 


0.5 


42 


(0,-1.5,-5.5) 


0(0.4) 


0.8 
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FIG. 6: Zero-field reduced integrated intensity ///T=80mK for 
tliree AF Bragg reflections as a function of reduced tempera- 
ture T/Tjv. The solid line is a guide to the eye. 



B. Magnetic order in applied field 



FIG. 7: Integrated intensity of the (0,0.5, -1.5) AF reflection 
vs. applied field at T = 80 mK. Intensity units are the same 
as in Table |. The solid line is a guide to the eye. The inset 
shows the Bragg peak intensity as a function of ^' (left) and 
2Q (right) in applied fields below (solid circles) and above He 
(open circles) . The solid lines correspond to a fit to the data 
in the (5', 2Q) plane as described in the text. 



The effect of an applied field on the magnetic order 
was investigated for fields along a. A spin-flop phase 
arises (spins cant out of the (&, c) plane towards the field 
axis) evidenced by a perpendicular antiferromagnetic or- 
der coexisting with a ferromagnetic moment along the 
field direction. Fig. shows the field dependence of the 
(0,0.5,-1.5) reflection measuring the antiferromagnetic 
component. The intensity increases with increasing field, 
reaches a maximum around ~ 1.4 T and then it 
abruptly drops to zero at He — 2.10(4) T in a sharp, 
near-first order phase transition. Throughout the spin- 
flop phase {0 < H < He) the ordering wave- vector was 
constant at the commensurate position (0, 1/2, 1/2). 

The transition at He was measured both with increas- 
ing and decreasing field, and no measurable hysteresis 
effect was observed as shown in Fig. 0. At 3 T, no 
AF or incommensurate reflections were observed along 
symmetry direction in the (6, c) plane with an inten- 
sity larger than 1.5% of the zero-field intensity of the 
strong (0,0.5,-1.5) reflection, equivalent to half the av- 
erage background level. The absence of magnetic Bragg 
peaks in the (6, c) plane suggests that the phase immedi- 
ately above the critical field He is a disordered phase (no 
LRO) stabilized by the applied magnetic fields c.f. the 
disordered SL phase in Fig. 1. 

Fig. H compares the field dependence of the reduced 
intensity I/Iq of three AF reflections at T = 80 mK, 
where Iq is the zero-field integrated intensity. The field 
dependence coincides showing that the magnetic struc- 
ture formed by the antiferromagnetic moments in the 
(6, c) plane is unchanged throughout the spin flop phase 
up to the critical field He- Furthermore, it shows that 
the initial increase in the AF Bragg peak intensities is 
not due to a rearrangement of the moments in a differ- 
ent spin configuration but arises from an increase in the 
magnitude of the antiferromagnetically ordered moment. 
This effect is attributed to the suppression of zero-point 
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IVlagnetic Field (TesIa) || a 

FIG. 8: Reduced integrated magnetic Bragg peak inten- 
sity, I/Io, as a function of field for three AF reflections at 
T = 80 K. The data for (0,1.5,1.5) and (0,0.5,-0.5) were 
measured with increasing field, while most of the data for 
(0,0.5, —1.5) was measured with decreasing field as shown in 
Fig. m. The solid line is a guide to the eye. 



quantum fluctuations by the applied field which allows 
more of the available spin moment to order. 

The effect of the applied field initially suppressing fluc- 
tuations and stabilizing the antiferromagnetic order is 
even more pronounced at elevated temperatures where 
the zero-field moment value is further reduced by ther- 
mal fluctuations. Fig. ^ shows that the antiferromagnetic 
moment increases significantly in applied field to reach a 
maximum around 1.4 T and only at much higher fields it 
collapses in a sharp transition at i?c=2.10(4) T (nearly 
T-independent up to 195 mK=0.9 T/v). 

Consistent with the above observation, scans in tem- 
perature at intermediate fields (below He) observe that 
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FIG. 9: Peak intensity of the (0,0.5, -1.5) AF reflection vs. 
applied field at three diff'erent temperatures. Intensity units 
are the same as in Fig. |l^. The solid lines are guides to the 
eye. 



CO 
CL 



CS2C0CI4 (0,0.5,-1 .5) H=1 .9T 



3000 
2500 
2000 
1500 
1000 



• T=85mK 
O T=250mK 




18.50 18.75 19.00 19.25 19.50 

2e 



100 150 200 250 

Temperature (mK) 

FIG. 10: Peak intensity of the (0, 0.5, -1.5) AF reflection vs. 
temperature in fixed applied field H = 1.9 T. Intensity units 
are the same as in Fig. H. The solid line is a guide to the eye. 
Inset: observed intensity of the (0,0.5,-1.5) reflection as a 
function of scattering angle 20 at two temperatures T = 80 
and 250 mK. The 250 mK data was scaled to match the peak 
intensity of the 80 mK data for direct comparison. The solid 
line is a Gaussian fit. 



the transition temperature is higher than in zero field, 
as expected for a structure with an increased ordered 
moment (stabilized by the field) that has a larger mean- 
field energy and is thus more stable against thermal fluc- 
tuations. This is shown in Fig. ^ by measurements of 
(0,0.5,-1.5) peak intensity in a field of T. The 

intensity decreases slower with increasing temperature 
than in zero field [see Fig. |^ and the order is stable 
beyond the zero-field transition temperature T/v=217(5) 
mK and disappears only at a much higher temperatures 
above the range covered by this experiment, extrapolated 
to Tn{H=1.^ T) ~ 300(20) mK. The long-range coher- 
ence of the structure is unchanged in the whole mea- 
sured temperature range as evidenced by the same angu- 



FIG. 11: Schematic magnetic phase diagram of CS2C0CI4 as 
a function of magnetic field and temperature. Solid symbols 
denote transitions where the AF LRO in the (6, c) plane disap- 
pears, see Figs. ^(|lo|) and ^(^ . Open symbols mark cross-over 
fields in the magnetization curve (see Fig. p^ ) identified with 
near-saturation of the lower doublet magnetization. Dashed 
lines are guide to the eye. 



lar width of 20 scans at the two extreme temperatures 
shown in the inset of Fig. lOl 



C. Temperature-field phase diagram 

In Fig. |l^ we show an {H \\ a,T) phase diagram of 
CS2C0CI4 based on the measurements described above. 
At zero field, CS2C0CI4 orders below rAr=217(5) mK 
with a commensurate wave- vector k = (0,1/2,1/2). 
The chains are ordered antiferromagnetically along their 
length and the moments are contained in the (&, c) plane 
at a small angle with the 6-axis. Magnetic fields applied 
along the a-axis initially stabilize the antiferromagnetic 
order by suppressing fluctuations. This is directly ob- 
served both in the increase of the perpendicular antifer- 
romagnetic moment in applied field (see Figs. ||-^) and 
also in the increase in the transition temperature at finite 
field (compare Figs. |3| and [lO|). At higher fields the order 
becomes unstable and above 7Jc=2.10(4) T|| a it collapses 
in a sharp, near-first order transition to a phase with no 
long-range magnetic order in the {b, c) plane, possibly a 
spin-liquid phase. Based on 1) the absence of order in the 
(&, c) plane, 2) nonsaturation of the magnetization, see 
Section HID, and 3) the expected order-disorder tran- 
sition driven by the large fluctuations arising from field 
noncommutation, we identify the phase above the criti- 
cal field He with the spin-liquid (S.L.) state predicted to 
occur below saturation in Fig. |l|. 



D. Ferromagnetic moment as a function of field 

The ferromagnetically-ordered moment was deter- 
mined from the intensity of the (Oil) reflection. The 
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FIG. 12: Ordered moments as a function of applied field 
along a at r=80 mK. Open circles show the perpendicular or- 
dered moment and solid symbols indicate the ferromagnetic 
(||a) component. The solid lines are guides to the eye and 
the dashed line shows the estimated partial contribution to 
the magnetization due to polarizing the higher-doublet states 
alone, assuming a constant susceptibility vs. field. The inset 
shows the canting angle 9 made by the total ordered moment 
with the (6, c) plane (above Hc=2.1 T it is assumed that the 
spin components perpendicular to the field axis are not or- 
dered and therefore 6—90°). 



resulting magnetization curve at T=80 mK is shown in 
Fig. |l^. The ferromagnetic moment increases over the 
whole range of the measurements up to 6.4 T and two 
regimes can be identified: a low- and a high- field re- 
gion separated by a cross-over around iJ„i=:2.52±0.06 T 
above which the rate of increase of the magnetization is 
significantly reduced. This cross-over is best illustrated 
in a plot of the differential susceptibility (x = dM/dH) 
in Fig. |l^(b) which shows significantly reduced values 
above Hm (of order 4±1 compared to zero field). This 
cross-over behavior is typical of spin-3/2 systems with 
two energe tically -separated Kramers doublets (for details 
see Section IV D ): in applied field the lower- lying doublet 
is saturated first at a field Hm above which the magne- 
tization curve shows a large decrease in susceptibility as 
only higher-doublet states can still be polarized. 

The cross-over field Hm (defined experimentally as the 
field where the magnetization is within less than 5% 
of the linear high-field behavior) is plotted in Fig. |ll| 
(open symbols). Note that the antiferromagnetic order 
in the (6, c) plane is suppressed at 7Jc=2.10(4) T much 
below the lower-doublet saturation field identified with 
iJ™=2.52±0.06 T (r=80 mK). The absence of long-range 
order for He < H < H„i is not due to thermal fluctua- 
tions as this field range remains finite after extrapolation 
to T=0 as shown in Fig. |ll|. This region is thus a quan- 
tum disordered phase induced by the applied field and 
is thus consistent with the proposed spin-liquid phase in 
the schematic phase diagram in Fig. |l|. 
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FIG. 13: (a) Magnetization vs. field along the a-axis at var- 
ious temperatures r=140, 180 and 250 mK. The solid line 
is a guide to the eye and the sloping dashed line shows a fit 
to a linear behavior at fields above ~3 T. (b) Susceptibility 
vs. field at r=80 and 140 mK. The solid line is guide to the 
eye (the horizontal dashed line shows the approximated par- 
tial contribution from the higher-doublet states). is the 
critical field where the antiferromagnetic order in the {b, c) 
plane disappears and Hm is the cross-over field above which 
the magnetization has approached the high-field near-linear 
behavior. 



IV. DISCUSSION 

In the previous section we have established the mag- 
netic structure of CS2C0CI4 and the effects of a noncom- 
muting field on its ground state. Here we discuss relevant 
microscopic mechanisms that may give rise to the ob- 
served ordered structure and find that a simple (nearest 
neighbor) mean-field picture can not explain all observed 
features. We then estimate the ordered moment and find 
a significant reduction from the available moment, indi- 
cating the presence of strong quantum fluctuations. Fi- 
nally, we relate the observed magnetic phase diagram to 
the minimal model of the XY-chain in a nonconimuting 
field. 



Commensurate vs. incommensurate order 



The observed ordering wave-vector condenses out of 
the diffuse scattering measured by Yoshizawa et alx\ in 
the disordered phase at 0.3 K just above T/v (inset of 
Fig. H), consistent with the expectation that order should 
arise at wavevectors where paramagnetie fluctuations are 
strong. The measured diffuse scattering3 showed system- 
atic incommensurate modulations along the chain direc- 
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tion attributed to competing interchain interactions. In 
isotropic systems, a mean-field picture for such frustrated 
couplings would predict incommensurate, spiral spin or- 
der along the chainsQ. The observed order however oc- 
curs at the commensurate, antiferromagnetic wavevector 
A;=(0,l/2,l/2) where interchain frustration effects cancel 
out and chains behave as decoupled. 

In CS2C0CI4 the tendency to form spiral order is sup- 
pressed because XY planes of neighboring chains are not 
parallel but make a large relative angle 2/3 resulting in 
an Ising-like frustrated interchain coupling between the 
XY-like ID chains. Considering only chains 1 and 4 in 
the chemical unit cell with XY spins the energy of a he- 
lical order in the easy planes of the two chains at a pitch 
to minimize frustration for the 6-axis spin components is 
Eheiix = S^[-J - |cos2/3|j2^/(2J)]. This spiral struc- 
ture becomes degenerate with the simple AF Neel or- 
der (spins along b) in the limit of orthogonal easy planes 
2/3=90° (close to the actual situation in CS2C0CI4). In 
the proximity of this limit of small effective frustration 
other effects may stabilize the observed commensurate 
Neel-type order and possibilities include: (1) zero-point 
quantum fluctuations could introduce non-linear terms 
in the free-energy expansion which may lift the classi- 
cal degeneracy and promote ordering at the decoupling 
point fc=(0,l/2,l/2), (2) further neighbor couplings such 
as between sites 1 and 5, or 1 and 6, although believed 
to be small could potentially stabilize the antiferromag- 
netic order, and (3) other effective Ising-type anisotropics 
arising from spin-orbit coupling or crystal field effects 
may favor spin ordering close to the 6-axis with the 
lowest energy achieved for an antiferromagnetic-type ar- 
rangement (constant moment on each site) as opposed to 
other incommensurate structures. We have also consid- 
ered dipolar couplings and found an increased energy by 
10^^ meV per spin for the observed order compared to 
a ferromagnetic-type arrangement along the c-axis and 
thus it was concluded that dipolar effects could not ex- 
plain the observed structure. 



B. Mean- field analysis 

To identify the exchange couplings involved in stabi- 
lizing the observed magnetic structure we calculate its 
energy in the mean-field approximation. We find sim- 
ple energetic arguments to explain: 1) why the observed 
magnetic structure belongs to the irreducible repre- 
sentation, 2) why the ordered moments are confined to 
the {b, c) plane and 3) why a multi-domain structure oc- 
curs. 

We first assume an isotropic Heisenberg exchange such 
that the three spin components along x, y and z (a- 
, b- and c-directions) could be treated separately. The 
observed ordering wave- vector fc=(0,l/2,l/2) indicates a 
doubling of the unit cell along the b and c axes leading 
to 16 different magnetic sublattices in the magnetic unit 
cell, as shown in Fig. 0(a). Using this extended unit cell 



we calculated the interaction matrix rj in the mean-field 
approximation; the eigenvectors of rj are basis vectors of 
the magnetic ordering and the eigenvalues are the cor- 
responding energy levels. Diagonalization of 77 gives the 
two lowest energy levels X± = —J±Jac, assuming weakly- 
coupled AF chains running along the 6-axis ( J > and 
J ^ Jab, Jac, Jbc)- For AF exchange between spins 1 and 
2 {Jac > 0) the ground state has energy A_ — — J — Jac 
and there are 3 degenerate eigenvectors: (1) ordering of 
the y{b) spin components in the irreducible representa- 
tion (Eq. pOl), (2) ordering of the z(c)-components 
also in F^'' or (3) ordering of the 2:(a)-components in F^ 
(Eq.|). 

Upon including anisotropy effects this three-fold de- 
generacy is lifted favoring ordering along the 5-axis {b is 
a common easy axis for all spins raS|-sliown by paramag- 
netic susceptibility measurementaldO) . The determined 
magnetic structure shown in Fig. ^ is indeed in the F^*^ 
representation with the largest spin moment along the 
easy 5-axis and a small moment along c, also in F^*'. The 
observed confinement of the ordered moments to the (6, c) 
plane can also be understood on energetic grounds: ac- 
cording to Eq. ^ ordering of spins along a in F^° would 
have parallel spins on sites 1 and 2, which would be en- 
ergetically unfavored by the AF interchain couplings Jac 
(the ordering of the a components can not belong to an- 
other representation, say F^, because the full Hamilto- 
nian including all anisotropy terms respects the symme- 
try of the crystal structure and therefore does not mix 
spin ordcrings from different irreducible representations). 

For a magnetic ordering of the b and c spin components 
in the F^*^ representation the spin configuration is not 
unique, but instead four distinct domains are possible, 
all with the same mean-field exchange energy. Those four 
domains are shown in Fig. ^ and the degeneracy arises 
because the F^" representation of the b and c components 
is two-dimensional (sites 1 and 3 and independent). 

Interestingly, at the mean-field level the ground-state 
energy does not depend on the couplings between sub- 
lattices (1,2) and (3,4). This is a general result for any 
antiferromagnetic ordering along b, i.e. fc=(0,l/2,Z) as 
can be easily seen by inspecting Fig. ^(a): each spin on 
the second group of sublattices interacts with pairs of 
antiparallel spins from the first group such that the net 
interactions cancel out. For example, spin 4 interacts 
through Jfcc with the pair of antiparallel spins 1 and 9 (5 
and 13). Similarly, spin 4 at 0.25 interacts through 
J^(, with the pair of antiparallel spins 6 and 14 at x^Q.lb 
{Jab with antiparallel spins type 6 and 14 at 0.25). 

The observed ordering at A:=(0,l/2,l/2) thus consists 
of two interpenetrating magnetic lattices that are non- 
interacting at the mean-field level. One lattice containing 
atoms 1, 2, 5, 6, 9, 10, 13, 14, and the other has atoms 
3, 4, 7, 8, 11, 12, 15, 16 in the magnetic unit cell. Since 
those two global lattices are related by inversion symme- 
try (1 at the center of the chemical unit cell) the intra- 
sublattice interactions are identical and so the ordered 
magnetic moment is expected to be the same for both, 
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as assumed in the analysis of the magnetic Bragg peaks 
in Section [II that gave good agreement with the exper- 
iment. 

The observed small alternating tilt = ±15(5)° of 
the magnetic moments away from the 6-axis indicates 
an effective local spin anisotropy. We do not have an 
explanation for its origin but propose that it may arise 
when the Heisenberg interchain exchange between sites 
with rotated XY planes (such as Jac between spins 1 and 
2 in Fig. |^) is projected onto the lower-lying Kramers 
doublet of effective spin-1/2. 



1. Full Hamiltonian 



C. Absolute magnitude of the ordered moment 



The magnitude of the ordered moment at T = 80 mK 
was determined by comparing the nuclear and the AF 
Bra gg p eak intensities (see Appendix C). Using Aexp and 
Eq. |ll| the resulting ordered moment at T = 80 mK is 
Too — 1.7(4)/i_B (zero field). 

The lower-lying Kramers doublet of Co^+ ions with ef- 
fective spin-1/2 has anisotropic g- values, gx,y = 2(7 and 
gz = g- Here 5=2.4 is the isotropic 5- value of the underly- 
ing spin-3/2 and was determined from high-tempjerature 
paramagnetic susceptibility measurementgli3ii3't3. This 
implies that the available moment in the XY plane is 
5a;,j/M_B/2=2.4 ^B- From the present diffraction experi- 
ments the estimated ordered moment along the in-plane 
axis b common to all spins is toq cos 0=1.6(4) fis^ clearly 
smaller than the available full moment in the XY plane 
(the ordered moment along the c-axis is very small and 
is a mixture between longitudinal and transverse parts). 
The reduction of the ordered moment from the full avail- 
able value indicates strong zero-point fluctuations in the 
ground state. 



D. Magnetization curve 



In this section we consider a minimal magnetic Hamil- 
tonian consistent with the crystal field and magnetiza- 
tion data on CS2C0CI4 and relate the observed magnetic 
phase diagram to the phenomenology of anisotropic mag- 
nets in noncommuting fields. We analyze the magnetiza- 
tion curve in terms of a spin-3/2 Hamiltonian appropriate 
for the Co^+ ions and identify the observed crossover at 
H„i with near-saturation of the lower- lying spin doublet. 
We then estimate the lower-doublet magnetization curve 
(Fig. |lj) and discuss it in terms of an effective 5=1/2 
XXZ Hamiltonian in noncommuting field. 



A minimal Hamiltonian for the spin 5=3/2 Co^+ ions 
that includes the ID exchange and crystal field effects is 



IS, ■ S,+, + D S:, 



- g^lBHJ^ - gfiBH^S^, 
(4) 

where / is the nearest-neighbor (isotropic) exchange in- 
teraction along the ID chains and D is the easy-plane 
anisotropy energy (perpendicular to the local z-axis). 
The last two terms in Eq. ^ are the Zeeman energy in lon- 
gitudinal {Hz = HsinP) and transverse {H^ = HsmP) 
fields, where (3 is the angle between the field direction 
and the XY plane (for H\\a this is equal to the angle 
between the local 2-axis and the c-axis in Fig. H). 

In the absence of magnetic fields {H = 0) and in the 
limit of large local anisotropy {D 3> /) only the lower- 
lying doublet I ± ^) contributes to the low-energy dynam- 
ics and in, this subspace the degrees of freedom can be 
describedEl by an effective spin 5=1/2 XXZ Hamiltonian 
(Eq. ^) with exchange J = 4/ and anisotropy parameter 
A=0.25. 



2. XXZ chain in a longitudinal (commuting) field 

The physics of the XXZ model (Eq. ^ is well under- 
stood in the limit A=0 when it is equivalent to a ID free 
fermion gas. Longitudinal fields {H \\ z) act as chemical 
potential filling up the quasiparticle band. The magneti- 
zation is directly related to, the filling factor and at finite 
temperature T is given bycl 



M{H,T) 



duj tanh 



gzfJ'sjH - Hiz cosa;) 
2fcBT 



(5) 

with the T=0 result M{H) = 2Mi^ sin"^ (iZ/iJi^,) /tt, 
H < Hiz. Fig. |lj shows a typical magnetization curve 
plot at T=Q (dashed curve) and finite T (solid line) where 
the approach to saturation is rounded off by thermal fluc- 
tuations. The saturation magnetization Miz = gzfJ.B/'2 
(per spin) is reached at the critical field (for finite A) 
Hiz = J{l + A)/gzfJ.B- The zero-field, T = longitudinal 
susceptibility is reduced compared to the semi-classical 
value due to i-auantum fiuctuations and for A=0.25 is 
calculated ast3 x z {0)^0 -237 g^ 11%/ J. Using the pro- 
posed values for the exchange interactions in CS2C0CI4, 
J=0.23(l) meV, gz=2A the longitudinal critical field is 
estimated as Hiz=2.1{l) T. 



3. XXZ chain in a transverse (noncommuting) field 

Transverse fields (i?||a;) have a very different effect 
compared to longitudinal fields because they : (1) break 
the spin rotational symmetry from U(l) to Ising, and (2) 
do not commute with the exchange terms. In low fields 
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this produces perpendicula]p,Jong-range antiferromag- 
netic order (spin-flop phase)Lnj. The non-commutation 
of the field introduces fluctuations that become strong 
enough above a critical field to suppress the antiferromag- 
netic order and induce a transition to a spin liquid state 
where the spin moment is not yet saturated along the 
field and the finite perpendicular spin components have 
exponentially-decaying correlations. A typical ground 
state phase diagram is shown in Fig. |^. The critical dis- 
ordering field is estimatedO at ~ J(3 -I- /^/ig^^B 
whereas exact diagonalizations of finite chainsO suggest 
that near-saturation of the magnetization occurs at much 
higher fields, similar to the classical saturation value 
Hix — 2J/gx^jLB- For CS2C0CI4 those estimates give 
~1.3(1) T and Hix^l.l{l) T. The zero-field, T = sus- 
ceptibility is again reduced compared to the semiclassical 
value due to quaptum fluctuations and for the XY chain 
is calculated afO Xxi'^) ^ 0.075glfi%/J. 



4- Full Hamiltonian in a longitudinal field 

For the S=3/2 Hamiltonian in Eq. H with well- 
separated energy scales for the lower | ± 2) and higher 
I ± |) doublets a cross-over in behavior is expected in ap- 
plied magnetic field from a low-field region where only the 
lower doublet participates (or both doublets depending 
on the field direction) to a higher field region where the 
lower doublet is saturated and only the higher-doublet 
contributes to the magnetization. 

Longitudinal fields {H \\ z) do not mix higher-doublet 
states until very large fields of the order the inter-doublet 
separation energy — 2D/ g^fJ^n- At low fields only 
the lower-doublet states contribute and the physics is 
that of the XXZ model in longitudinal fields. Above the 
lower-doublet saturation at Hiz the magnetization shows 
a plateaux at Mi^ = gz^J-s/^ (per spin) stable up to fields 
around when the magnetization increases again by 
mixing in states from the higher-lying doublet | ± |) to 
finally reach the full spin value of 3/2giJ.B- The upper 
longitudinal critical field for CS2C0CI4 is estimated at 
^2^=9(1) T. 



5. Full Hamiltonian in a transverse field 

Transverse fields {H \\ x) in Eq. ^ have finite matrix 
elements between the two doublets and thus mix higher- 
doublet states into the ground state (of order g^sH /2D) 
at any finite field.E^I Both doublets participate at low 
fields (large susceptibility) and above a cross-over field 
the lower doublet is near saturated and only the higher 
doublet contributes (small susceptibility). Such a behait. 
ior is evident in earlier magnetization measurementSLj'E£l 
on CS2C0CI4 in fields along the 6-axis (entirely trans- 
verse) where the magnetization increases rapidly at low 
fields with a large susceptibility then crosses over above 



T to a region where the magnetization increases much 
slower to approach saturation (3.6 hb) at fields >16 T. 



6. Comparison with the observed (total) magnetization 

The a-axis magnetization shown in Fig. ^ is in broad 
agreement with the expected behavior for mixed longi- 
tudinal and transverse fields on the full S=3/2 Hamilto- 
nian in Eq. ^: a rapid increase is observed at small fields 
(both doublets contribute) followed by a cross-over to a 
much slower increase at higher fields (only higher-doublet 
states contribute). The cross-over field i/„i=2.52±0.06 
T (r=:80 mK) is identified with near-saturation of the 
lower doublet. This field appears to be larger than the 
estimated near- or full-saturation fields for purely trans- 
verse Hix or purely longitudinal fields Hiz, possibly due 
to either the approximations used in estimating Hi^ solv- 
ing Eq. ^ assuming decoupled lower and higher doublets, 
or other terms in the Hamiltonian, such as interchain 
exchanges not included explicitly here. 

Fig. lU shows the partial lower-doublet magnetization 
assuming the contribution to the magnetization from the 
higher doublet states can be approximated by a con- 
stant susceptibility (as shown in Fig. [l^ ). This assump- 
tion is consistent with calculations in the single-site ap- 
proximation following Ref. |l^, which predict that the 
higher-doublet susceptibility for fields tilted at (3=45° is 
relatively small (compared to Xx,z{Q) of the lower dou- 
blet) and decreases smoothly by only ~25% between zero 
and 5 T. Putting the magnetization values in Fig. |lj 
on an absolute scale gives a lower-doublet saturation 
moment Ms=1.7(4) ^b, consistent with typical values 
expected for fields applied at an angle to the z-axis 

Ms = '^g^cos^P + g^ sin^ /3/is/2 ~1.9 hb for /3=45°. 



7. (Partial) magnetization of the XXZ chain 

Since no detailed predictions are available for the mag- 
netization curve of XXZ chains in mixed longitudinal and 
transverse fields we compare the results with the generic 
analytic form given in Eq. ^. This is valid strictly only for 
XY chains in longitudinal fields, but it provides a simple 
analytic form to parametrize the data and extract an ef- 
fective saturation field , the saturation magnetization 
and an effective "temperature" T, which is a measure of 
the fluctuations causing the rounded approach to satura- 
tion. 

The extracted saturation field i?5=2.37(3) T (80 
mK<r< 250 mK) is signiflcantly larger than the criti- 
cal field He— 2.10(4:) T where the antiferromagnetic order 
disappears, giving further support for the existence of an 
intermediate phase between the antiferromagnetic order 
and the nearly fully-polarized phase as indicated in the 
phase diagram in Fig. |l|. The fitted "temperatures" f 
are systematically larger than the true measuring tem- 
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FIG. 14: Partial magnetization from lower-doublet states 
alone as a function of applied field along a. The data is ob- 
tained from the observed total magnetization by subtracting 
an estimate of the higher-doublet contribution as indicated in 
Fig. ^ (dashed line) and described in the text. Hc=2.1 T 
is the critical field where the antiferromagnetic order in the 
{b, c) plane is suppressed (see Fig.^ and Hm is the cross-over 
field above which the (lower-doublet) magnetization is nearly 
(within 5%) saturated. The solid line is a fit to the generic 
form in Eq. n (the dashed line shows this calculation at zero 
temperature) as described in the text. 



peratures T (T=200(50) niK for the r=80 niK data, fit 
shown by solid line in Fig. indicating more fluctua- 
tions in the system than can be accounted for by tem- 
perature alone. A source of those fluctuations can be the 
noncommuting fields that create disorder effects at large 
fields. At small noncommuting fields the dominant effect 
is breaking the spin rotational symmetry which has the 
consequence of promoting long-range order in a spin-flop 
phase. This is also in agreement with the experiments, 
which observed that low fields stabilize the perpendicular 
antiferromagnetic order. 



E. Further studies 

The measurements presented here have highlighted a 
potentially very significant field induced phase transition 
in CS2C0CI4. The field dependence of the order corre- 
lates well with the expectations of a quantum magnet 
driven through a quantum critical point by a noncom- 
muting field. Essential to this picture is the identification 
of the disordered phase with a gapped spin liquid state. 
Further studies using NMR, magnetic susceptibility, and 
heat capacity are called for to investigate whether a spin 
gap really does open above the critical field of 2.10 T. It 
would also be very interesting to look for any evidence of 
spin glass behavior induced by the field. Inelastic neu- 
tron scattering measurements of the excitations in a field 
should be particularly revealing and we plan to make 



such measurements in the near future. 



V. CONCLUSIONS 

In conclusion, single-crystal neutron diffraction was 
used to determine the magnetic ordering as a function of 
noncommuting applied field in the quasi-lD spin-1/2 XY- 
likc antifcrromagnet CS2C0CI4. In zero field long-range 
order with wavevector fc=(0,l/2,l/2) was found below 
T/v=217(5) mK. The magnetic structure was determined 
using group theory and has spins ordered antiferromag- 
netically along the chains with moments confined to the 
(6, c) plane. A domain structure was found with adjacent 
chains in different phases. Possible mechanisms promot- 
ing this commensurate order were discussed. The ob- 
served reduction in the ordered moment was attributed 
to zero-point fluctuations in the ground state. 

Magnetic fields applied along the a-axis were found to 
initially stabilize the perpendicular AF order and form a 
spin-flop phase. This structure becomes unstable at high 
fields where a transition occurs to a phase with no LRO 
in the (5, c) plane. Measurements of the ferromagnetic 
component found that near-saturation of the moments 
occurs only at much higher fields. The phase in-between 
the spin- flop and saturated phases 2.10 < Hsl < 2.52 T 
II a has been proposed to be a spin liquid state disordered 
by the strong quantum fluctuations arising from the ap- 
plied field noncommutability with the exchange Hamil- 
tonian. 
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Appendix 

A. Group theory analysis 

Using the space group symmetry of the crystal struc- 
ture we identify allowed basis vectors for a magnetic 
structure with the observed wavevector fc=(0,l/2,l/2). 
This is done by determining the irreducible representa- 
tions and eigenvectors of the little group Gk of symmetry 
operations that leave the wavevector k invariant. 
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We start by considering the symmetry elements of the 
Pnma space group of CS2C0CI4: 

{1, T, }, (6) 

where 1 is the identity operator, 1 is the inversion at 
the origin, 2q, denotes a 180° screw axis along direction 
a = x,y or z (180° rotation followed by a translation 
with half a unit cell along axis a) and where through- 
out this group-theory section the axes x, y and z refer 
to the crystallographic directions a, h and c. niap is a 
glide plane containing axes a and (3. The group is non- 
symmorphic because the group elements {i?|a} consist of 
an operation R followed by a translation a equal to half 
a direct lattice vector. The experimentally observed or- 
dering wave- vector k = (0, 1/2, 1/2) is invariant under all 
these operations such that the little group Gk contains 
all of the above elements. 

The representations of the little group Gk are given by 
exp(— ifca) r"(i?fe). The ordering in CS2C0CI4 is a spe- 
cial case because the crystal symmetry is nonsymmorphic 
and the ordering wave-vector lies on the Brillouin zone 
boundary. The representations of elements Rk of the 
space group with subsequent translation t = (0,n,TO) 
such that n + m = odd have a different sign to those 
of operations Rk alone, and this leads to additional ir- 
reducible representations. They can be found, by adding 
new group elements to the little group Gfcjl3 which are 
the original elements with an additional translation t. We 
call the additional translations 

{1*, T*, 2l 2l, 2*, mly, mL, J . (7) 

Thus the little group Gk consists of a total of 16 ele- 
ments. We determined the classes and the character table 
of this group and these are shown in Table ||. The group 
consists of 10 different classes and therefore has 10 irre- 
ducible representations. Only two of the representations 
fulfill the necessary condition that x(l) = ~x(l ) (this 
condition follows directly from the prefactor of the rep- 
resentation which changes sign under a translation by t). 
The experimentally observed ordered magnetic structure 
is thus associated with one of these two representations. 
The eigenvectors (j/^ of the irreducible represeptations 
were determined using the projector method.ES They 
are given by 

= A = 1...10, (8) 

g 

where g is an element of the little group and (p is any 
vector of the order parameter space. For fc=(0,0.5,0.5) 
the magnetic unit cell is doubled along the h and c axes, 
and the order parameter space is a 48-dimensional axial 
vector because the magnetic unit cell contains 16 mag- 
netic moments as indicated in Fig. ||)(a) and each has 
three space components. We obtain 

</>^ = (TOla;,'Tllj^,mi2,-TOl:r,mij/,mi2, 



-mi, -7712, -mz, -rrii, 
-nil, -7712, -ms, -mi, 

mi, 1712, ms, 1124) (9) 

for representation F^. mia is the component a of the 
magnetic moment i with 1 < i < 16. For better 
readability, in the above equation the individual com- 
ponents of the magnetic moments 5 to 16 were omitted 
and the magnetic moments were written as vectors, i.e. 
ms — —mi = {—mix,—miy,—miz). For representation 
r^'^ we obtain 

0io = {mix, miy, mi^, mi^, -miy, -mi^, 
"^3a; , rri'iy , m^^ , m^^ , -msy , -msz , 
-mi, -m2, -m-i, -m^, 
-mi, -m,2, -m-i, -m^, 

mi,m2,m-i,mi) . (10) 

The ordered magnetic structure is thus doubly degener- 
ate in each of the three spin components given by the 
dimension 2 of these two representations. 

A magnetic structure in the irreducible representa- 
tion can occur in several different domains. Assume spins 
are confined to the (6, c) plane in a typical configuration 
shown in Fig.^(a) called domain Al. One can construct 
domain A2 in Fig.^c) by reversing the 6(?/)-components 
of all spins and keeping the c(z)-components unchanged. 
Domain Bl in Fig|^(b) is obtained from domain Al by 
reversing the spins on sites 3, 4 (and 7, 8, 11, 12, 15, 
16) and using the same rule one transforms domain A2 
in Fig.|(c) into B2 in Fig.|(d). Domains Al and A2 
are indistinguishable from each other in a neutron scat- 
tering experiment because they have the same magnetic 
structure factor and domains Bl and B2 are also indis- 
tinguishable. However, an A-type domain (either Al or 
A2) has a different structure factor from a B-type (either 
Bl or B2) domain. 



B. Zero-field magnetic structure 

The integrated intensity of magnetic Bragg peaks is 
related to the structure factor of the magnetic ordering 
througftJ 

,^Q,.{p.)\J^^\lm,,M)?. (11) 

\2fiBj Kio sm(2e) 

where f{Q) is the magnetic form factor for Co^"*" ion 
$ is the flux of incident neutrons, and Vmo are 
the number and the volume of the magnetic unit cells, 
7 = 1.913 and tq = 2.818 • 10"^^ m. I{Q) is the to- 
tal integrated intensity of a Bragg reflection measured in 
the (^',20) plane and l/sin(20) is the Lorentz correc- 
tion factor that arises because intensity is measured as 
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TABLE II: Irreducible representation of the group Gk- 



a function of angular coordinates. F±(Q) is the compo- 
nent of the magnetic structure factor perpendicular to 
the scattering wave vector and is defined as 



FAQ) = FiQ)-{F{Q)-Q)Q, 

where Q is the normalized wave- vector transfer. 
The magnetic structure factor is defined as 



(12) 



HQ) 



16 

E 



nii exp(-i Q ■ di) 



(13) 



where Q is the wave- vector transfer in the experiment, 
di are the positions of the Co^^-ions and the sum is over 
all 16 magnetic ions in the magnetic unit cell. For a two- 
domain structure (A -|- B type) the magnetic Bragg peak 
intensity can be written as 



I{Q)^alA{Q) + {l-a)lB{Q). 



(14) 



where a is the population of the A-domain and (1 — a) 
is the population of the B-domain. 



C. Absolute magnitude of the ordered moment 

The magnitude of the ordered moment at T = 80 mK 
was determined by comparing the nuclear and the AF 



Bragg peak intensities. The intensity of a nuclear Bragg 
peak is given as 



Vb sin(2e) 



(15) 



where TV and Vq are the number and the volume of the 
unit celk, Fm{Q) is the nuclear structure factor and 
given asll3 



Fn{Q) = ^6iexp(-iQ • di) 



(16) 



where the sum runs over all elements in the nuclear unit 
cell and bi is the elastic scattering lengthEd of atom i in 
the unit cell. The nuclear peaks used for calibration were 
(022), (033), (022) and (033) and their observed relative 
integrated intensities were consistent with the calculated 
structure factors to within 15%. Multiple scattering and 
extinction corrections were assumed to be negligible. The 
measured nuclear intensities gave the overall scale factor 



for the intensities 



N 



(2^) 



$ in Eqna. The mag- 



nitude of the ordered magnetic moment was determined 
using ^oxp E^nd Eq. |ll|. The resulting ordered moment at 
T — 80 mK in zero field is too — 1.7(4)/iB. 
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